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Abstract 

We present an analysis of the cosmological evolution of matter 
sources with small anisotropic pressures. This includes electric and 
magnetic fields, collisionless relativistic particles, gravitons, antisym¬ 
metric axion fields in low-energy string cosmologies, spatial curva¬ 
ture anisotropies, and stresses arising from simple topological defects. 
We calculate their evolution during the radiation and dust eras of 
an almost isotropic universe. In many interesting cases the evolu¬ 
tion displays a special critical behaviour created by the non-linear 
evolution of the pressure and expansion anisotropies. The isotropy 
of the microwave background is used to place strong limits of order 
< 5 X 10“®A (1 -|- Zrec)~^ on the possible contribution of these 
matter sources to the total density of the universe, where 1 < A < 3 
characterises the anisotropic stress. The present abundance of an 
anisotropic stress which becomes non-relativistic at a characteristic 
low-energy scale is also calculated. We explain why the limits obtained 
from primordial nucleosynthesis are generally weaker than those im¬ 
posed by the microwave background isotropy. The effect of inflation 
on these stresses is also calculated. 
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1 Introduction 


One of the key problems of cosmology is deciding which matter helds are 
present in the Universe. The most notable uncertainties govern the presence 
or absence of influential matter helds which might contribute to the total 
density of the Universe today, and thereby offer a resolution of the ’dark 
matter’[1] problem or provide an effective cosmological constant or vacuum 
stress [2] which alters the signihcance of the evidence governing the process 
of large-scale structure formation. Other matter helds, like magnetic helds 
or topological defects, can have a profound inhuence upon the evolution and 
properties of galaxies. In the early universe the identity of matter helds is 
even more uncertain. String theories are populated by large numbers of high- 
mass helds but they will not survive to inhuence the late-time evolution of 
the Universe signihcantly if inhation takes place at high energies. Inhation, 
of course, involves it own material uncertainties: it requires the existence of 
a weakly-coupled scalar held which evolves slowly during the early universe 
[3]. In all these cases the principal problems are ones of fundamental physics 
- do the postulated matter helds exists or not?, what are the strengths of 
their couplings to themselves and to other helds? Relic densities residing in 
isotropic or scalar stresses are the simplest to determine. They can be cal¬ 
culated from a semi-analytic solution of the Boltzmann equation if their in¬ 
teraction strengths, lifetimes and masses are known. Limits on their masses, 
lifetimes, and interaction cross-sections are obtained by requiring that they 
do not contribute signihcantly more than the closure density to the Universe 
today, or produce too many decay photons in a particular portion of the 
electromagnetic spectrum [1], [4]. 

In this paper we are going to consider the behaviour of anisotropic stresses, 
of which electric and magnetic helds are the simplest examples. We shall not 
assume that inhation has taken place, although we will also calculate the 
ehects of inhation on their present abundances. The isotropy of the mi¬ 
crowave background requires any such anisotropies to be extremely small 
but we shall hud that their evolution is subtle because of the coupling be¬ 
tween pressure anisotropies and the expansion dynamics during the radiation 
era. Anisotropic stresses have many possible sources — besides cosmologi¬ 
cal magnetic and electric helds [5], examples are provided by populations of 
collisionless particles like gravitons [6], photons [7], or relativistic neutrinos 
[8], [9]; long-wavelength gravitational waves [10], [11], Yang-Mills helds [12], 
axion helds in low-energy string theory [13], [14]; and topological defects like 
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cosmic strings and domain walls [15]-[17]. Our experience with high-energy 
physics theories also alerts us to the possibility that a (more) final theory 
of high-energy physics will contain many other matter fields, some of which 
may well exert anisotropic stresses in the universe. 

Cosmological observations of the isotropy of the microwave background 
tell us that the Universe is expanding almost isotropically. Therefore we can 
expect to employ the high isotropy of the microwave background to constrain 
the possible density of any relic anisotropic matter fields in the Universe. We 
shall see that a very general analysis is possible, largely independent of the 
specific identity of the physical field in question, which leads to strong limits 
on the existence of homogeneous matter fields with anisotropic pressures, 
irrespective of initial conditions in many cases. This analysis exhibits a 
number of interesting features of the general behaviour of small deviations 
from isotropy in expanding universes in the presence and absence of a period 
of inflation. 

In section 2 we set up the cosmological evolution equations for a universe 
containing a perfect fluid and a general form of anisotropic stress. Under 
the assumption that the anisotropy is small these equations can be solved. 
They reveal two distinctive forms of evolution in which the abundance of 
the anisotropic fluid is determined linearly or non-linearly, respectively. The 
evolution can be parametrised in terms of a single parameter which charac¬ 
terises the pressure anisotropy of the anisotropic fluid. In section 3 we give a 
number of specific examples of anisotropic fluids and explain how this analy¬ 
sis also allows us to understand the evolution of the most general anisotropic 
universes with 3-curvature anisotropy with respect to the simplest examples 
which possess isotropic 3-curvature. In section 4 we derive limits on the 
abundances of anisotropic stresses in the Universe today by using the COBE 
four-year anisotropy data [18]. We also consider the case of an anisotropic 
stress which becomes non-relativistic after the end of the radiation era (anal¬ 
ogous to a light neutrino species becoming effectively massive). In section 5 
we consider the effect of a period of de Sitter or power-law inflation on the 
abundances of anisotropic stresses surviving the early universe, and in section 
6 we compare the constraints we have obtained on anisotropic stresses from 
COBE with those that can be derived by considering their effects on the pri¬ 
mordial nucleosynthesis of helium-4. The effect of the pressure anisotropy is 
to slow the decay of shear anisotropies to such an extent (logarithmic decay 
in time during the radiation era) that in the general case the COBE limits 
are far stronger than those provided by nucleosynthesis. 
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2 The Cosmological Evolution of Skew Stresses 

Consider an anisotropic universe with metric [19] 


ds"^ = dt^ — a?{t)dx^ — h‘^{t)dy^ — (?{t)dz^. 
If we define the expansion rates by 

^ a ^ a ^ 
a = -; P = 7; ^ 

a 0 c 

where overdot denotes d/dt, then the Einstein equations are 


( 1 ) 


cx -\- 0 ? oii^jd T — 

-s^g{t; - ir), 

(2) 

(3 -\- (3^ -\- (3(oi + ^) = 

-87rG(T| - It), 

(3) 

e + + e{a + p) = 

-87rG(T| - ir), 

(4) 

a b c 

“ + I + “ = 

a b c 

-87rG(T“ - ir), 

(5) 


where is the energy-momentum tensor and henceforth we set SttG = 1. 
We define the mean Hubble expansion rate, H, by 

H = ( 6 ) 

and the two relative shear anisotropy parameters by 

R = —^ and ^ ^ . (7) 

When i? = S' = 0 the universe will be the isotropic flat Friedmann universe. 

We are interested in studying the behaviour of a non-interacting combi¬ 
nation of a perfect fluid and a fluid with an anisotropic pressure distribution. 
Thus the energy-momentum tensor for the universe is a sum of two parts: 

K ^tt + si ( 8 ) 

where is the energy-momentum tensor of the perfect fluid, so 
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= dmg{p„ -p., -p„ -p,) 


( 9 ) 


with equation of state 

p* = (7 - l)p*, (10) 

and is the energy-momentum tensor of a fluid with density p and principal 
pressures Pi,P 2 , and ps, so 

si = diag{p, -pi, -p 2 , -ps). (11) 

We shall assume that the principal pressures of the anisotropic fluid are each 
proportional to its density, p, so that we have 

(Pi,P2,P3) = (Ap,i^p,/ip), (12) 

where A, u, and p are constants. It will be useful to define the sum of the 
pressure parameters by 


/S. = \ + p + v. (13) 

In the special case of an axisymmetric metric, (0)> two of these three 
constants would be equal. The ratio of the densities of the two fluids is 
defined by 

C = -. (14) 

P* 

The two stresses tl and si are assumed to be separately conserved, so 


the a 
tions. 


VX = 0 


V^s' 


b’) 


( 15 ) 


0 components of which, using (p!0D-(p!^, give two 


conservation equa- 


P* + 3hf7p* = 0, (16) 

p -|- p[a(A -|- 1) -h /3(/i -h 1) -|- Oiy -\- 1)] = 0. (17) 

Hence, the perfect fluid density falls in proportion to the comoving volume, 
as 




1 

{abc)^ ’ 


( 18 ) 
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and the anisotropic flnid density falls as 

One can see from these two eqnations that Q oc admits 

a variety of possible time-evolntions according to the valnes of the isotropic 
and anisotropic pressnres. The evolntion eqnations become, 


cn -j~ cn -jr cr(/3 T — 

/3 + + (5{oi + 0) = 

e + e'^ + e{a +13) = 


-^[P + Pi - P2 - P3 + 

1. 

^[P + P3-Pl-P2 + 

1. 

-[P + P 2 -P 3 -P 1 + 


(2 - 7)P*] 

(20) 

do 

1 

* 

(21) 

(2 - 7)P*] 

(22) 


Snbstitnting H,R and S for a, (3, and 9, from (|^) and nsing (^0|) -(^ ) and 
(H), we have two propagation eqnations for the anisotropies, 


HR + RH + 3RH^ = Pi — Ps = {X — ^)Qp^, (23) 

HS + SH + 3SH^ = Pi - P2 = (A - z/)gp„ (24) 

and the conservation eqnations (H)-® combine to give 


Q = ^|9(7 


1) — 3(A + u + p) + R{2fi — n — A) + S{2v — \ — p)] (25) 


The three eqnations (|2^)-(|25|) completely determine the evolntion when the 
anisotropy is small. We are interested in the solntions of these eqnations in 
the case where the anisotropy is realistically small. Thus we assume that H 
and take the values they would have in an isotropic Friedmann universe 
containing isotropic density p^ (ie. S' = i? = Q = 0), so 


2 4 

H = -and p* = — 

3-it 3'yH^ 


(26) 


This is consistent with the ([]) — equation, (^. The three equations (p3|)- 
then completely determine the evolution of R, S, and Q. They reduce to 
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( 27 ) 


R+-A2-',) 

it 

s + -i2-i) 

it 

Q 


2Q(A-/i) 

it 

2Q{X - p) 
it 

20 

-^^[ 9(7 — 1 ) — 3A + R{ 2 fi — u — X) + 
+S{ 2 i' — A — /i)]. 


(28) 

(29) 


Axisymmetric solutions exist when S' = 0 and X = v. No essential simplifi¬ 
cation arises by imposing axial symmetry and so we shall treat the general 
case with S' 7 ^ 0. 

The system of equations 0-0 gives rise to a characteristic pattern of 
cosmological evolution when stresses with anisotropic pressures are present. 
Equation (|5I1|) reveals that there is a critical condition which, if satished, 
makes the problem a second-order stability problem; that is, if we linearised 
the equations about the isotropic solution (i? = S' = 0) with zero anisotropic 
stress density {Q = 0) we would hnd a zero eigenvalue associated with the 
evolution of Q{t). This critical condition is 

Criticality condition : 3(7 — 1 ) = A. (30) 

When 37 < 3 -|- A the shear distortion variables R and S relax towards their 
attractor where i? = S = Oast—^-cx), and so we have 


R = 


S = 


2Q{X-y) 

2-7 


2-7 


■ + 8it T- 

, 2 g(A-/i) 

(31) 


2-7 

+ 82 t '•< - 

, 2Q(A-Z/) 

(32) 


2-7 

7-2 




with (5i and 82 constants. The 8 it~ terms are the contributions from the 
isotropic part of the 3-curvature. They fall off more rapidly than the part of 
the shear that is driven by the anisotropic pressure (assuming A 7 ^ /i 7 ^ z/). 
As t grows, the 8 terms become negligible if 37 > 2A (and recall that the 
isotropy is stable so long as 37 < 3-|-A). Thus, R and S become proportional 
to Q(t), which is determined by 
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Q^^{9(7-1) 


■3A + 


2Q 


(2-7) 


[(A -/i)(2/i - 1 /- A) + (A - z/)(2z/- A -/i)]}. 

(33) 

When the evolution is not critical ( 3(7 — 1) 7 ^ A) the evolution of Q{t) 
is determined by the terms linear in Q on the right-hand side of (|33D; eqns. 
(lUD’dll) sfill hold, but now we have 


Q{t) = Qot^; (34) 

k = ^[ 3 ( 7 - 1 )-A]. (35) 

37 

Our assumption, eq.(^), that the evolution of and H follow their values 
in an isotropic Friedmann universe will only be consistent at large times if 
k < 0, that is if 


-3 < 3 ( 7 -1) < A < 3. (36) 

The right-hand side of this inequality derives from the causality conditions for 
signal propagation in the i = 1, 2, 3 directions (pi < p so A < 1, p < 1, z/ < 1); 
the left-hand side arises from pi > —p, which ensures the stability of the 
vacuum. When k > 0 the anisotropic stress redshifts away more slowly 
that the isotropic perfect fluid on the average (over directions) and comes 
to dominate the expansion dynamics, making them completely anisotropic. 
We do not live in such a universe. By contrast, when k <0, the isotropic 
stresses redshift away the slowest and increasingly dominate the dynamics, 
so the gravitational effect of the anisotropic stresses steadily diminishes. In 
the critical case the average stress energy of the anisotropic stress is counter¬ 
balanced by the isotropising effect of the perfect fluid and the shear evolution 
is determined by the second-order effects of the pressure anisotropy. This 
effect can be seen in the study of free neutrinos by Doroshkevich, Zeldovich 
and Novikov [8] and in the study of axisymmetric magnetic fields by Zeldovich 
| 20 ], 

When the evolution is critical, (|30|) holds, k = 0, and the evolution of 
Q{t) is decided at second-order in Q by 




4g" 

97(2 — 7 )t 


{(A — p)(2p — z/ — A) -|- (A — z/)(2z/ — A — p)}. (37) 






Hence, 


Q{t) 

A 


_Qo_ 

1 + AQo In 


(38) 


- /^)(2/i - z/ - A) + (A - iy)(2iy - A - /i)}, (39) 


where Qo and ti are constants. For physically realistic stresses we have H > 0, 
and so as t —>• cx) the ratio of the energy densities approaches the attractor 

while the associated shear distortions approach the values given by eqns. 

). Since the values of R and S at the epoch of last scattering of 
the microwave background radiation determine the observed temperature 
anisotropy we will be able to constrain the allowed value of Q at the present 
time by placing bounds on R and S at last scattering and then evolving the 
bounds forward to the present day. 

We can also calculate the asymptotic forms for the expansion scale factors 
of (^. In the critical cases, as f —>• cx), they evolve towards 




a{t) 

oc t^jlnt}’". 

(41) 

b{t) 

oc t^{lnt}"'. 

(42) 

c{t) 

oc t^jlnt}”. 

(43) 

where the constants m, n, w are 

dehned by 


m = 

4(2A — jji — v) 

97 / 1 ( 2 - 7 ) ’ 

(44) 

n = 

4(2/i — V — \) 

(45) 

97 / 1(2 - 7 ) 
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w 


4(2z/ — n — X) 
97^(2 — 7 ) 


(46) 


Some specific cases will be considered below. Note that the asymptotic form 
for the scale factors, (^TD - (^3D is consistent with the principal approximation 
imposed by eq. (|^) . 

In the non-critical cases (with A; < 0), as t cx), the scale factors evolve 
to hrst order as 


a{t) 

(X - Wt*’}, 

(47) 

b{t) 

oc 

(48) 

c{t) 

oc t^{l-Vct^}, 

(49) 


where the constants I 4 , 14 , and I 4 are dehned by, 

{ 14 , H, Vc} = ^ {h + ^ ~ 2 A, a + /i — 2iy, X + — 2/i}. (50) 

These asymptotes reveal the different behavionr in the critical cases which 
prodnces the logarithmic decay of the shear. In the non-critical cases the 
anisotropic pertnrbations to the scale factors decay as power laws in time 
since A: < 0. 


3 Some Particular Skew Stresses 

There are many examples of anisotropic stresses to which the analysis of the 
last section might be applied. We shall consider some of the most interesting. 
In each case we can determine the characteristics of the ’critical’ state in 
which the evolntion of the anisotropy is determined by the non-linear conpling 
with the anisotropic stress. In section 4 we will go on to calcnlate the observed 
microwave backgronnd temperatnre anisotropy. The most familiar example 
of an anisotropic stress is that of an electromagnetic held. 
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3.1 Magnetic or Electric fields 

In this case the anisotropic energy-momentum tensor of eqn. ( 0 ) has a 
simple form [20], [22], [24]. For example, for a pure magnetic held of strength 
B directed along the z-axis, we have Tq = = —Tl = —T 2 = i?^/ 87 r and 

so this corresponds to the choice 

X = u = = 1 (51) 


Hence, 


^{magnetic) = A{electric) = 1 


(62) 


Therefore, from eqn.(pOD, we see that the criticality condition is obeyed 
when the background universe is radiation dominated. Thus, in the standard 
model of the early universe, the evolution of magnetic (or electric) helds 
will exhibit the non-linear logarithmic decay found in equations (|T7D-(|40D. 


This was hrst pointed out by Zeldovich [20], and can be identihed in the 
calculations of Shikin [21] and Collins [22]. Barrow, Ferreira, and Silk [23] 
used the COBE data set [18] to place constraints on the allowed strength of 
any cosmological magnetic held. 

As a specihc example of a critical case, a radiation dominated universe 
(7 = 4/3) containing a magnetic held aligned along the z-axis, X = v = 
—fi = 1, gives A = 4 in (|39D and hence 


a{t) (X 6 (f) cx {lnt }4 ( 53 ) 

c(f) cx f 2 {lnf }“2 (54) 

Notice that the volume expansion goes as in the isotropic radiation universe, 
abc cx f 2 in accord with the principal approximation stipulated by (|26|) . 

As an example of a non-critical case consider a pure magnetic held aligned 
along the ; 2 -axis of a dust universe (7 = 1). We have k = —2/3, and so at 
late times 


a{t) 

c{t) 


cx 6(f)cxf3{l-^ 


3f 3 


cx 




3fi 
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This behaviour can be seen in the exact magnetic dust solutions of Thorne 
[24] and Doroshkevich [25]. Other studies of the evolution of cosmological 
magnetic helds in anisotropic universes can be found in refs. [28]. 

The axion held in low-energy string theory [13], [14] also creates stresses 
of this characteristic form and a source-free magnetic held has been used to 
study the possibility of dimensional reduction in cosmologies with additional 
spatial dimensions near the Planck time by Linde and Zelnikov [26] and 
Yearsley and Barrow [27]. 

3.2 General Trace-free stresses 

The magnetic held case is just one of a whole class of skew helds which exhibit 
non-linear evolution during the radiation era. Any anisotropic stress which 
has a trace-free energy-momentum tensor, s^, will have 

A = 1 (55) 

and will exhibit critical evolution in the presence of isotropic black body radi¬ 
ation with 7 = 4/3. This case includes free-streaming gravitons produced at 
tpi ~ 10 “^^s which are collisionless at f because of the weakness of the 

gravitational interaction mediating graviton scatterings [ 6 ]. It is also likely 
to include all asymptotically-free particles at energies exceeding lO^^GeV 
when interparticle scatterings, decays and inverse decays have interaction 
rates slower than the expansion rate of the universe, H. In all cases dehned 
by A = 1 the shear to Hubble expansion rate decays only logarithmically 
during the radiation era. 


R (X S (X 



(56) 


In the dust era the critical condition will not continue to be met for trace- 
free fields. Although the evolution of Q, R, and S is now determined at linear 
order in (|^), and (|40|), there is still a significant slowing of the decay of 
isotropy by the anisotropic stresses compared to the case where anisotropic 
stresses are absent. We see that when 7 = 1 we have 


— X R X S X 
P* 


1 

W 


(57) 
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compared to a decay of i? oc S' oc when anisotropic stresses are ab¬ 
sent {Q = 0) or ignored. The study by Maartens et al [29] of the evolution 
of anisotropies in a dust-dominated universe containing possible anisotropic 
stresses at second order (described by a tensor in reference [28] which is 
equivalent to the used here) displays this same slowing of the shear decay 
(the discussion of ref. [41] omits this consideration). 

In general, we can see that if the criticality condition is satished in the 
radiation era, when 7 = 4/3, then it cannot be satished during the dust era 
that follows, when 7 = 1 . However, an interesting situation can arise if the 
source of the trace-free stress is a population of particles which are relativistic 
above some energy (so A = 1 there) and non-relativistic when the universe 
cools below this energy scale (so A = 0 there). Thus there can be a change 
in the value of A with time. We shall examine this case in section 4. 

Some matter helds with anisotropic pressures, like Yang-Mills helds [12], 
correspond to an anisotropic huid with time-dependent A, /i, and u. But 
if the evolution is slow enough, it is well-approximated by the model with 
constant values of A, /i, and z/ used here. 

3.3 Long-wavelength gravitational waves 

In the last section we considered the evolution of anisotropic stresses in the 
simplest flat (Hq = 1) anisotropic cosmological model of Bianchi type I with 
the metric ([^). The most general anisotropic universes of this curvature 
are of Bianchi type VIIq and they have an Einstein tensor that can be 
decomposed into a sum of two pieces: one corresponding to the Einstein 
tensor for the simple type I geometry, the other to a piece that describes 
additional long-wavelength gravitational waves [10], [11], [30]. These waves 
create anisotropies in the 3-curvature of the universe in addition to the sim¬ 
ple expansion-rate anisotropies present in the Bianchi I universe (which has 
isotropic 3-curvature). However, the contribution by the long-wavelength 
gravitational waves can be moved to the other side of the Einstein equa¬ 
tions and reinterpreted as an additional ’effective’ energy-momentum tensor 
describing a ’fluid’ of gravitational waves. It has vanishing trace. This de¬ 
composition means that any general flat (or open) Bianchi type universe of 
type VIIq (or Vllh) containing an isotropic perfect fluid, ( 0 ) behaves like 
a Bianchi type I (or type V) universe containing that fluid plus an addi¬ 
tional traceless anisotropic fluid. The parameters A, /i, and u are approxi¬ 
mately constant when the anisotropies are small. Hence the general evolution 
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of anisotropic universes with anisotropic curvature containing isotropic ra¬ 
diation will exhibit the same characteristic logarithmic decay of the shear 
anisotropy given by eqns. ®. ( 0 . and (^OD during the radiation era. 
This behaviour appears in Collins and Hawking [31], Doroshkevich et al [32], 
and other authors, [33], [34], During the dust era any anisotropic curvature 
modes will behave like trace-free stresses and, although the evolution will no 
longer be critical, the shear anisotropy will fall more slowly {R oc S oc ) 
than in simple isotropic universes with isotropic curvature and no anisotropic 
stresses (where R oc S oc 

3.4 Strings and Walls 

Topological defects provide another class of anisotropic energy-momentum 
tensors which fall into the category of stress modelled by (|^ for part of their 
evolution. The energy-momentum tensors of string sources were also consid¬ 
ered more generally by Stachel [16], Marder and Israel [17]. The specihc de¬ 
scription of line stresses created by topological defects is reviewed in ref. [15]. 
An inhnite string with mass per unit length fi extending in the x-direction 
contributes an anisotropic stress-tensor = ^S{z)6{y)diag{l,1,0,0)] that 
is 


A = —1, /X = z/ = 0 — A{string) = —1. (58) 

Their evolution could therefore only be critical in a universe containing 
a perfect fluid with equation of state p* = —p*/3. It is worth noting that 
this corresponds to the evolution of the curvature term in the Friedmann 
equation when the universe is open. Thus we would expect the string stress 
to evolve critically during the late curvature-dominated stage {1 + z < Dq 
of an open universe, during a curvature-dominated pre-inflat ionary phase, or 
during any period when quantum matter helds with 7 = 1/3 dominate the 
expansion. 

For slow-moving inhnite planar domain walls of constant surface density 
in the x — y plane , the stress corresponds oc g{z)diag{l,1,1,0), where 
g{z) is a local bell-shaped curve centered on 2 ; = 0 [15],and corresponds to 
the choice 


A = n = —1, p = 0 — A{waU) = —2. 


(59) 
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The wall stress would only be critical if the equation of state of the back¬ 
ground universe is p* = — 2p^/3. 

These descriptions do not include the complicated effects of the non¬ 
linear evolution of a population of open strings and loops which must include 
intersections, kinks, gravitational collapse and gravitational radiation. 

3.5 General classification and energy conditions 

When evaluating the form of the residual anisotropy and energy density 
created by anisotropic cosmological stresses it is most convenient to classify 
stresses by the value of A. We can circumscribe the likely range of realistic 
A values by considering some general restrictions on the energy-momentum 
tensor. The dominant energy conditions [35] require us to impose the physical 
limits 


|A| < 1, |/i| < 1, |z/| < 1 . (60) 

Hence we have the bounds 

-3<A<3. (61) 

If the strong-energy condition [35] were imposed on the energy-momentum 
tensor then we would have a stronger restriction 

A > -1. (62) 

This condition is violated by string and wall stresses and necessarily by any 
isotropic stress which drives inflation (0 < 7 < 2/3) since A{wall) = —2 
and A{string) = —1. We see from (|34[)-(^^ that there can only be approach 

to isotropy at late times in the radiation era if A > 1. When 0 < A < 1 

the expansion approaches isotropy during the dust era but not during the 
radiation era. When A < 0 isotropy is unstable during both the dust and 
radiation eras. 


4 Microwave background limits 

Let us consider the simplest realistic case in which the Universe is radiation 
dominated until l + Zeq = 2.4 x then dust dominated thereafter; 

here; hlo A 1 is the cosmological density parameter and Hq is the present 
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value of the Hubble constant in units of Mpc~^. We shall assume 

that the microwave background was last scattered at a redshift Zreo where, 
in the absence of reheating of the cosmic medium, 

1 + Zrec = 1100. (63) 

If there is reionization of the universe then last scattering can be delayed 
until 1 + Zrec = 39(Qhoho)~^ for Qo^rec << 1, where Qho is the present baryon 
density parameter [36]. 

The microwave background temperature anisotropy is determined by the 
values of shear distortions R and S at the redshift Zrec- The evolution of the 
photon temperature in the x, y, and 2 ; directions is given by 


T 

T 

T 


Tn 


Tn 


Tn 


Oo 

a{t) 

bp 

b{t) 

Cq 

c{t) 


Toexp{—y adt}, 

(64) 

Toexp{- J pdt], 

(65) 

Toexp{— J 6dt}. 

(66) 


If we dehne the temperature anisotropy by 


iT (T, - T,) + (% - T.) 

T To 


( 67 ) 


then, for small anisotropies (so exp{—/ adt} ~ 1 — / adt} etc), using the 
dehnitions of (0) in (|67D , we have 


— = -Hj(R + S)dt (68) 

Since recombination always occurs in the dust era, H = 2/3t,and the ob¬ 
served microwave background anisotropy will be 


x/('d,0,Ho)- (69) 

where / ~ 0(1) is a pattern factor taking into account non-gaussian statis¬ 
tical factors and the possible pattern structure created by complicated forms 
of anisotropy in the general case [13], [23],[37]. Typically, in the most gen¬ 
eral homogeneous flat universes the pattern combines a distorted quadrupole 
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with a spiral geodesic motion. In ref. [23] a detailed discussion of the sam¬ 
pling statistics of the microwave background temperature distribution and 
the non-gaussian nature of the anisotropic pattern was given. This leads to 
bounds on / in flat and open universes of 

0.6 </ < 2.2. (70) 

The discussion above shows that the A > 1 case is the realistic one which 
allows evolution towards isotropy at late times. The observed anisotropy is 
therefore given in terms of the present value of the density ratio, Q{to), by 

6T 2Q{to)f{2X-fi-u){l+Zrec)^ 

T A ^ ^ 

If we take the COBE 4-year data set to impose a limit of 6T/T < 10“® on 
contributions by the anisotropic fluid, and use ([70|) , then the present density 
parameter of the anisotropic fluid, hlao, is limited by 


f^aO A 


8.3 X 10-®A 


A 


(2A - /i - Z/)(l -f Zrec)' 


2A — /i — z/ 


1100 
1 T 


X 10 


-5.08-3.04A 


(72) 


Since 1 < A < 3 for realistic fluids we can examine the extremes of 
this limit which is strongest when there is no reionization of the Universe at 
2 ; << 1100 because reionization allows a longer period of power-law decay 
of R, S, and Q, hence a smaller residual effect on the microwave background 
isotropy. In the two extreme cases we have limits of 


A = 1 : 12,0 < 3.7 X 10-^ X ( ™ ) (73) 

\ 1 “T ^rec ^ 

/ linn \3 

A = 3 : < 9.5 X 10“^^ X -) (74) 

\ 1 ~h 

The A = 1 case with the choice of ( 0 ) corresponds to limits on the present 
cosmological energy density allowed in magnetic (or electric helds) which was 
studied in [20] and in [23] where the most general evolution of anisotropy was 
included. Note that these limits are far stronger than those that are generally 
obtained for isotropic forms of dark matter by imposition of astronomical lim¬ 
its on the maximum total density, the age of the Universe, or the deceleration 
parameter [1], [4]. 
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4.1 Evolution of anisotropic dark matter with a char¬ 
acteristic energy scale 

We are familiar with the standard picture of the cosmological evolution of 
light (m << IMeV) weakly interacting particles [1], [4]. When T > m 
they behave like massless particles with their number density similar to that 
in photons, up to statistical weight factors. Their energy density redshifts 
away like {1 + zY until the temperature falls to the value of their rest mass. 
Thereafter they behave like non-relativistic massive particles and their en¬ 
ergy density redshifts away more slowly, as (1 -|- We can consider an 
analogous scenario for anisotropic stresses. Suppose we have an anisotropic 
fluid which behaves relativistically until the temperature falls to some value 
T+ < ~ leV and then behaves non-relativistically at lower tempera¬ 

tures. This situation corresponds to following the evolution with A = 1 for 
T > T+ and then with A = 0 for T < T+. The evolution is complicated 
by the fact that it will be critical during the radiation era but non-critical 
during the first stage of the dust era when the temperature exceeds T+ when 
Q will decay in accord with (^), before becoming critical again during the 
remainder of the dust era until the present when Q will decay logarithmi¬ 
cally in accord with (^Ol). In this case the final density of the anisotropic 
dark matter is constrained by the microwave background isotropy to have 


VL, 


\ao < 1.2/-' X 10-^ 

\ 1 “t“ Z^^f^ J 


where 


1 

^ Aln(l Zj^) 


(75) 


A = -{(/u — A)(2/i — ly — A) + (ly — A)(2z/ — A — /u)} (76) 

9 

For example, if the anisotropic matter has a characteristic energy scale m 
then, since m = = 2.4x 10-^(l-l-z+) eV, we have (picking A = n = — /i = 1 


so A = 32/9), with (|70D , that, 

flao < 5.7 X 10-5 (^: 


1100 


AAZrec) (8.3 -|- InMev) ^ ^ 

and there is a very weak dependence on the mass scale m^v = (m/lel7). For 


m > leV we have, roughly, that 

Qao < 6.9 X 10-5 


1100 

1 ~h ^re 


(78) 
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Therefore these fields are always constrained by the microwave back- 
gronnd anisotropy to be a negligible contribntor to the total density of the 
Universe. 


5 The Effects of Inflation 

So far we have ignored the conseqnences of any period of inflation in the very 
early nniverse [3]. The eqnations (^7[) - (|3TD hold for the case of generalised 
(power-law) inflation with 0 < 7 < 2/3, [37]. However, for de Sitter inflation 
with 7 = 0 , the isotropic density drives the inflation with = 3Hq = 
constant, and they are changed to 

R + 3RHq = 3HqQ(\ — p), (79) 

S + 3 SH 0 = 3HoQ{X-iy), (80) 

Q = -^^[-9-3A + i?(2p-i/-A) + 5(2i/-A-p)].(81) 
3 

This system can only be critical for the evolntion of Q only if A = —3 
bnt this cannot occnr for an anisotropic flnid snbject to (|BD|). As t —*• cxo we 
have 

R —> q{X — p) + 5i exp{—3Hot) (82) 

S q{X — u) + S2exp{—3Hot) (83) 

Q (5oexp[-(3A)i7o^] (84) 

Thns, since A -|- 3 > 0 the contribntion made to the shear by the anisotropic 
pressure decays. Moreover, we see that the contribution by the isotropic 
curvature (5) terms to the shear dominates the contribution by the pressure 
anisotropy at late times if A > 0. In both cases the shear anisotropy decays 
away exponentially fast. This is in accord with the expectations of the cos¬ 
mic no hair theorems [38] because the anisotropic stresses considered here 
obey the strong energy condition when A > —3. Thus if N e-folds of de 
Sitter inflation occur in the very early universe, then the values of the shear 
distortion parameters, R and S, are each depleted by a factor exp(—3A) 
whilst the relative density in the anisotropic fluid, o, is reduced by a factor 
exp{-(3 +A)A}. 
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Similarly, in the case of power-law inflation (0 < 7 < 2/3), the isotropic 
curvature mode will dominate the late time evolution of the shear during the 
inflationary phase if 


7 < 


2A 


(85) 


If anisotropic stresses can be generated at the end of inflation then the 
analysis of the previous sections applies. 


6 Primordial Nucleosynthesis 


It is instructive to consider the question of whether primordial nucleosyn¬ 
thesis can provide stronger limits than microwave background isotropy on 
the densities of anisotropic stresses in the Universe. The key issue that de¬ 
cides this is whether the microwave background constraint on the expansion 
anisotropy created by an anisotropic energy density, which is of order 10 “® 
and imposed at a redshift Zrec ~ 10 ^ or lower (in the event of reionization), 
is stronger than a limit of order 1 — 10 “^on changes to the mean expansion 
rate imposed at the epoch of neutron-proton freeze-out, Zfr r'U 10^°. The is¬ 
sue is decided by knowing how fast the anisotropic stresses decay with time 
between Zfr and Zrec- 

In the simplest anisotropic universes, which have isotropic 3-curvature, 
the shear anisotropy falls like the that of the 5 mode in eqns. (|^ - (|3^ . 


Consider hrst the simple isotropic curvature case with no anisotropic stress, 
so Q = 0. The evolution of the shear to Hubble rate parameters, R and S', 
during the dust and radiation era of a universe with flo = 1 is 


> Zeq '■ R <x. S (X t (X. 1 + Z ( 86 ) 

2 ; < Zeq '■ R X S X t~^ CX (1 -|- z)^ (87) 

Hence, if nucleosynthesis gives an upper limit of ~ 0.2 (roughly equiv¬ 
alent to adding one neutrino type to the standard three-neutrino model) 
on the values of |i?| and [S'! at the redshift of neutron-proton freeze-out, 
Zfr, this corresponds to an upper limit at the time of last scattering of 
0.2(1 -I- Zfr)~^{l + Zeq) ^(1 -|- Zrec)^, SO we have 
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rsj 


3 

2 


( 88 ) 


Since the microwave background gives a limit of 6T/T ~ (i? + S)rec < 
10“®,[18], we see that the nucleosynthesis limit is much stronger when the 
3-curvature is isotropic, as hrst pointed out by Barrow [39]; nor are conceiv¬ 
able improvements in microwave receiver sensitivity ever likely to close to gap 
of 10^ that exists between the nucleosynthesis and microwave limits on the 
isotropic 3-curvature 6 modes. These effects of these simple anisotropy modes 
on nucleosynthesis were considered in the papers of Hawking and Tayler [40] 
and Thorne [24], 

However, the situation changes when we consider the most general forms 
of anisotropy with anisotropic curvature or when matter is present with 
anisotropic pressures. If the evolution will be critical during the radiation 
era the anisotropy falls off so slowly during the period before Zeq that the 
microwave background limits become stronger than the nucleosynthesis lim¬ 
its. This effect was considered in both the dust and radiation eras by Barrow 
[39] in the evolution of Bianchi type VII universes close to isotropy but the 
observational limits were 100 times weaker then. For example, in the in¬ 
teresting case where A = 1 the evolution of anisotropies follows the form. 


•( 13 ), 


z > Zeq-. Roc S oc (89) 

< Zeq '■ R OC S OC 1 + Z. (90) 

Hence, the nucleosynthesis limit translates into an upper limit on R and S 
at the redshift of last scattering of only 

Rrec ~ ^r-ec < 7.1 X 10-"(Hoh2)-' X • (91) 

We see that this is never stronger than the microwave background limit of 
< 10“® for any possible redshift of last scattering. Alternatively, we might 
restate this result as follows: it is possible for anisotropic fluids to create a 
measurable temperature anisotropy in the microwave background radiation 
without having any signihcant effect upon the primordial nucleosynthesis of 
hehum-4. 


21 






7 Discussion 


We have set up the cosmological evolution equations for a very general class 
of anisotropic stresses in the presence of an isotropic perfect fluid. Such 
anisotropic stresses encompass a wide range of physically interesting cases, 
including those of cosmological electric and magnetic helds, a variety of 
topological defects, gravitons, and other populations of collisionless parti¬ 
cles. When the nniverse is close to isotropy the mean expansion rate behaves 
as in the isotropic Friedmann nniverse to leading order, bnt the density of 
the anisotropic stresses is conpled to the expansion anisotropy in an inter¬ 
esting way. It was fonnd that there are two possible forms for the evolution. 
In the critical case the density of anisotropic matter is determined by its 
non-linear coupling to the expansion anisotropy and density of anisotropic 
stress falls only logarithmically relative to the isotropic backgronnd density. 
By contrast, in the non-critical case the density falls as a power-law in time 
relative to the background density. In all cases the microwave backgronnd 
anisotropy can be used to place a limit on the density of matter that could 
be residing in the nniverse today in forms with anisotropic pressures. These 
limits arise because their evolution approaches an asymptotic attractor in 
which the temperature anisotropy produced in the microwave backgronnd 
by the anisotropic stresses is simply related to their density relative to the 
isotropic background density. The limits obtained on the possible cosmolog¬ 
ical density of matter of this form are far stronger than conventional limits 
on isotropic forms of dark matter because they make use of the microwave 
isotropy limits {ST/T < 10“^) rather than the far weaker limits from obser¬ 
vations of the Hubble flow and age of the Universe (Hq < We went on 

to explain why the limits that can be impose upon anisotropic stresses by the 
microwave backgronnd isotropy measnrements are in general stronger than 
those arising from the limits on the primordial nucleosynthesis of helium-4. 
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